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ABSTRACT 

Linear  and  nonlinear  studies  are  made  of  two-dimensional  free-surface 
flows  under  gravity,  in  which  a  disturbance  is  caused  to  an  otherwise-uniform 
stream  by  a  distribution  of  pressure  over  the  free  surface.  In  general,  such 
a  disturbance  creates  a  system  of  trailing  waves.  However,  there  are  special 
disturbances  that  do  not,  and  some  categories  of  such  disturbances  are 
discussed  here.  This  work  has  potential  applications  to  design  of  splash-less 


ship  bows. 
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SIGNIFICANCE  AND  EXPLANATION 


One  of  the  important  problems  of  modem  ship  hydrodynamics  concerns  the 
flow  at  the  extreme  bow.  In  particular,  reduction  and,  if  possible, 
elimination  of  the  splash,  and  resulting  splash  drag  is  of  great  interest. 

In  the  present  paper  linear  and  nonlinear  studies  are  made  of  two- 
dimensional  flows,  in  which  an  otherwise  uniform  stream  is  disturbed  by  a 
distribution  of  pressure  over  the  free  surface.  It  is  found  that  some 
families  of  pressure  distributions  do  not  generate  waves.  These  wave-less 
flows  generate  candidate  shapes  for  splash-less  bows. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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1 .Introduction 


Ono  of  the  important  problems  of  modern  ship  hydrodynamics 


concerns  the  flow  at  the  extreme  bow.  In  particular. 


prediction  and.  If  possible,  elimination  of  the  splash,  and 


resulting  splash  drag  component  is  of  great  Interest.  It 


appears  (c.f.  refs  (1,2))  that  one  of  the  roles  of  a  bulbous 


bow  Is  to  effect  Just  such  an  elimination  (or  at  least 


reduction  In  magnitude)  of  the  splash  drag  under  suitable 


conditions  of  speed  and  loading. 


In  previous  research  (3,4]  the  present  authors  have 


considered  fundamental  questions  associated  with  a 


two-dimensional  model  for  such  bow  flows.  That  work  suggested 


that,  for  a  given  bow  shape,  splash-free  solutions  may  not 


exist.  However,  only  special  (specifically  non-bulbous)  bow 


shapes  were  used  In  the  previous  studies.  The  possibility 


exists  that,  by  considering  suitable  famlJles  of  bow  shapes. 


one  could  Identify  a  special  shape  having  the  splash-free 


property.  In  the  linearised  case,  such  shapes  have  been 


considered  by  Schmidt  [3]. 
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In  the  prtunt  irtlclt,  we  approach  thla  problta  in  an 
Indirect  Banner.  If  ouch  a  apodal  bow  exlata,  and  If  ono 
cooput o •  tho  prosouro  exortad  on  It  by  tho  water,  than  ono  can 
replace  the  original  body  by  an  artificial  preecrlbed  preaeure 
dlatrlbut Ion  over  the  f rec-aur f ace .  Converaely,  If  one  eterte 
with  a  given  f ree-aur f ac e  preaaure  dl a t r lbut Ion ,  .the 
atreaellnea  generated  will  define  a  candidate  bow  ahape. 

The  obvloua  difficulty  with  thla  approach  la  that.  In 
general,  auch  preaaure  diatrlbutlona  generate  a  train  of  wavea 
downstream,  and  we  are  (ultimately)  looking  for  flowa  that  are 
not  only  aplaah-free,  but  alao  wave-free.  However,  it  la  well 
known  (c.f.{6],  p.404)  that  in  two  dlaenelona  there  exlat 
apeclal  preaaure  diatrlbutlona  whoae  downatreaa  wave  amplitude 
la  xero,  and  our  aaln  purpoae  here  la  to  enumerate  and  diacuae 
auch  apeclal  caaea. 

In  particular,  auppoae  we  conalder  a  family  of  preaaure 
diatrlbutlona  that  vanlah  for  x  >  0,  and  tend  to  a  poeltlve 
conatant  value  aa  x  -  -«  •  Thla  preaaure  will  in  general 
correapond  to  a  aem  1-inf  in  1 1  e  "item"  profile  having  finite 
draft  at  x  ■  and  producing  a  train  of  wavea  downatreaa  at 
x  ■  Within  thla  family  there  may  be  one  or  more  apeclal 
membera  for  which  tha  downatreaa  wave  amplitude  la  zero.  Since 
there  are  no  wavea,  the  radiation  condition  doea  not  apply  to 
auch  apeclal  aolutlona,  and  the  flow  direction  can  be  reveraed, 
ao  generating  a  aeml -inf Ini te  'bow"  flow.  Thla  flow  poaaeaaea 
neither  a  aplaah  nor  (upatraaa)  wavea,  and  hence  la  precleely 
of  the  type  dealred. 


The  above  discussion  proceeds  without  any  conalderatlon  of 
linearisation.  Hovavar,  faw  demonstrations  hava  baan  made  of 
non-linear  fraa-ourface  flowa  without  downatreaa  waves,  and  It 
la  not  entirely  obvious  that  they  exist,  although  the  numerical 
evidence  (7,8]  la  strong,  and  the  approximate  non-llnear  theory 
of  Tulin  (9)  can  be  Interpreted  aa  providing  Indirect  analytic 
confirmation. 

Here  we  first  discuss  the  linearised  case  In  some  detail, 
enumerating  some  families  of  wave-less  pressure  distributions, 
and  computing  the  corresponding  free-surface  contour.  This 
contour  can  then  be  replaced  by  a  fixed  body,  but  only  as  an 
approximation  when  that  body  has  a  small  stream-wise  slope 
everywhere,  in  order  that  the  generating  pressure  can  be  small 
to  justify  linearisation.  Some  of  the  results  suggest  bulbous 
character ,. but  of  course  a  true  bulb  can  never  be  attained  with 
such  a  restriction  on  the  body's  slope.  These  linear  results 
can  be  viewed  as  an  Inverse  re-formulation  of  the  theory  of 
Schmidt  (3]. 

We  then  solve  the  non-llnear  problem  by  a  perturbation 
procedure.  In  which  the  linearised  solution  Is  viewed  as  the 
first  term  In  an  asymptotic  expansion  with  respect  to  some 
small  parameter  e.  Thus,  a  measures  either  (and  both)  the 
body's  streamwlse  slope  and  the  size  of  the  pressure 
distribution  that  generates  it.  In  principle,  providing  we 
adopt  a  careful  Inverse  formulation  with  the  velocity  potential 
as  independent  variable,  any  number  of  terms  in  the  asymptotic 
expansion  In  powers  of  e  can  be  computed,  for  continuous 
pressure  distributions.  We  provide  here  some  2nd-order 


numerical  rtauld,  and  higher-order  raaulca  to  be  not  hard 
to  obtain.  Thla  la,  howavar,  not  trua  If  tha  Input  praaaure 
poaaaaaaa  step-f unction  dlacontlnultlaa ,  bacauaa  a lngul ar i 1 1  a a 
occur  vhoaa  degree  lncraaaaa  with  tha  powar  of  a. 

An  altarnatlve  approach  to  non-llnaar  frea-aurfaca 
aolutlon  la  via  axact  lnvaraa  aathoda,  c.f.  (10,11).  That  la. 
It  la  poaalbla  to  wrlta  down  a  apaclf  lcatlon  for  tha  flow  In 
which  all  requisite  aquations  (Including  tha  non-llnaar 

fraa-aurface  condition)  ara  satisfied  exactly.  Thla  haa  not  in 
tha  paat  proved  a  vary  uaaful  technique  bacauaa  of  tha 
unraallatlc  nature  of  tha  atraanllnaa  generated.  Howavar,  tha 
extra  freedou  obtained  by  allowing  a  non-taro  fraa-aurface 
praaaure,  enablea  ua  to  get  a  little  cloaar  to  practical 
applicability .  Tha  price  paid  for  thla  la  that  wa  cannot 
expect  tha  lnveraely-obtalned  praaaure  dlatr lbutlone  to  be 
exactly  taro  anywhere  on  tha  free  surface.  Tha  beat  wa  can  do 
la  to  causa  tha  pressure  to  be  vary  snail  over  a  portion  of  the 
atrearullne  that  la  a  candidate  for  a  trua  "free"  surface. 
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2.  Week  Pressure  Distributions 


Suppose  Chat  Che  only  disturbance  Co  a  uniform  stream  of 
unit  magnitude  In  the  ♦x-dir«ction  consists  of  a  small 
departure  P(x)  from  the  equilibrium  atmospheric  value  of  the 
free-surface  pressure.  If  P(x)  is  sufficiently  small,  and  we 
make  the  usual  assumptions  [12]  of  classical  linearised 
water-wave  theory  we  find  (c.f.  [13])  that  the  free-surface 

elevation  la  y  ■  n(x),  where 

nCx)  -  -P(  — >  «-J  P  *  ( t)  K  ( x-f  )  dt.  (2.1) 

The  kernel  function  K(x)  is  the  wave  elevation  due  to  a  unit- 
step  pressure,  namely 


K(x) 


^f(x)  ♦ 


'  2cos  x— 1 ,  x  >  0 

0  ,  x  <  0, 


(2.2) 


where 


f(x)  “  Cl(|x|)sln  x  -  sl(|x|)cos  x  sgn  x  (2.3) 


le  the  odd  auxiliary  function  ([14],  p.232)  of  the  sine  and 

cosine  Integrals. 


Note  that  we  have  non-dimensional  1  sed  so  that  not  only  Is 
the  free  stream  of  unit  magnitude,  but  the  acceleration  of 
gravity  and  the  density  are  both  also  of  unit  magnitude.  In 
effect,  we  heve  chosen  a  length  scale  of  magnitude  U*/g,  and  a 
pressure  scale  of  magnitude  pU*,  where  U,  g  and  p  are  the 
actual  values  of  these  quantities.  Since  the  wavelength  of 
linearised  waves  Is  2sU*/g,  this  means  that  in  the  present 


m 


non -d l*«m t on* 1  franework,  their  length  la  2a,  aa  la  clear  from 
the  tera  In  coa  x  In  (2.2). 


In  fact,  alnce  f(x)  tenda  to  zero  at  Infinity,  aa  x  - 
we  have 


n(x)  -  -P(  — )  + 


P  '  (  £  )  l  2  coa  (x-£)-l)d£ 


-  -P(+»)  +2  P' (£ )coa(x-()d( . 


(2.4) 


However,  froa  now  on  we  aaauae  (without. loaa  of  generality) 
that  P(+«)  ■  0,  i.e.  that  the  preaaure  la  allowed  to  return  to 
ataoapherlc  far  dovnmtream  %  and  hence  that  the  lean  free  — 
aurface  level  downatreaa  la  defined  to  be  the  level  at  which 
y  ■  0.  In  that  caae,  (2.4)  reducea  to 


-  "I' 


P'(()coa  £  d£ ]  coa  x  + 


P'(£)aln  £  d£ J  aln  x 


(2.5) 


l.e.  there  la  a  wave  at  x  ■  whoaa  aagnitude  and  phsee  are 
deteralned  by  the  preaaure  dletrlbutlon  P(x).  Thla  wave 
dlaappeara  If  and  only  if 


P'(x)coa  x  dx  ■  0 


(x)  aln  x  dx  •  0. 


(2.6) 


(2.7) 


Any  preaaure  dlatlbrutlon  P(x)  aatlafylng  (2.6)  and  (2.7)  **  ** 
be  deacrlbed  aa  "wave-leea" . 
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Per  exaaple,  auppeif  P ( x )  Is  step-wise  constant,  l.a. 

P<x)  “  P  j  "  constant,  Xj.,  <  x  <  Xj  (2.*) 

for  soae  aet  of  points  xj,  J»0 , 1 , 2 ,  .  .  .  ,  S ,  where  we  aay  take 
xe  ”  — •  and  xo  ■  0.  Than  the  preaaure  gradient  la  a 
collection  of  Dirac  S-functlona,  l.e. 

N 

P’(x)  -  I  (Pj+,-Pj)a(x-xp  (2.9) 

J- ‘ 

where  P(j+,  *  0.  Now  (2.6)  and  (2.7)  require 
N  N 

I  (Pj*.,-  P j ) coe  xj  •  t  ( P j  +  t  -  P j)aln  xj  -  0.  (2.10) 

J-i  j-» 

We  can  fix  one  of  the  conatanta  P j  i  aay  P,  ;  there  are  then 

(2N-2)  free  parametera  P2  ,  P . .  and  x,,x2 . , 

and  Just  two  equatlona  restricting  thea.  Clearly  there  la  no 
wave-lass  preaaure  If  N  ■  1.  Thla  alaply  conflraa  that  a 
(single)  step-function  pressure  always  generates  a  trailing 
wave  of  non-sero  amplitude. 

If  N  ■  2,  there  are  as  aany  equatlona  as  unknowns  and  a 

wave-less  pressure  Is  possible.  The  specific  fora  taken  by 
(2.10)  is 

(P,-P,)«o.  *,-?!-  ( P  »-P  t ) . in  x,  -  0  (2.11) 

Thus,  xt  “  -s,  -2s,  -3s,...,  so  that  the  step  length  must  be 
an  integer  aultlple  of  s.  If  that  Integer  Is  odd,  then  P,  ■  0 
and  P|  ■  Pt/2,  l.e.  the  pressure  P(x)  consists  of  a  juap  at 

x  •  x,  to  exact ty  hmlf  of  its  (non-xero)  upatreaa  value, 

followed  by  a  Juap  to  xero  at  the  origin.  In  practice,  we  are 


r^V^.  A  V 


■•Inly  interested  in  the  shortest  step,  with  x,  ■  -s,  whi.h 
corresponds  to  •  Froude  nusber  (based  on  the  step  length)  of 
*-«/*  >  .56*. 


On  the  other  hand,  the  solutions  with  x,  •  -2», 
-*s,  -6s,...  aust  have  P,  ■  0,  i.e.  correspond  to  a  pressure 
distribution  that  vanishes  outside  the  finite  segment 
x,  <  x  <  0  and  is  constant  within  that  segment.  The  existence 
of  such  staple  finite  wave-less  pressures  is  well  known;  for 
example,  the  aost  important  case  xt  -  -2s,  with  a  Froude 
number  of  .399,  was  discussed  by  Lamb  (6)  and  used  by  Schwarts 


The  absence  of  waves  in  this  case  is  easy  to  understand 
intuitively,  since  there  is  exactly  one  wavelength  (2s) 
between  beginning  and  end  of  the  pressure  distribution.  Since 
these  ends  generate  exactly  equal  and  opposite  waves, 
destructive  interference  is  possible.  In  fact,  a  similar 
argument  also  explains  the  semi-inf  lnl  to  case  x,  »  -s,  where 
the  step  length  is  exactly  half  a  wavelength.  Now,  in  view  of 
the  fact  that  the  pressure  steps  down  by  a  factor  of  1/2  at 
both  ends,  the  waves  produced  by  the  ends  are  exactly  equal  in 
magnitude  and  sign,  but  now  are  shifted  by  a  half-wavelength, 
so  that  again  destructive  interference  occurs. 


In  summary,  pressure  distributions  consisting  of  Just  two 
distinct  steps  can  generate  wave-less  free  surfaces,  providing 
the  distance  between  the  steps  is  an  integer  multiple  of  the 


half-wavelengi 


The  "shortest"  such  wave-less  pressure  has 


a  step-length  of  t ,  i.e.  one  half-wavelength,  and  corresponds 


-\V . V--  -VA  -\  -\ 


to  a  ■•■i-inf lnlte  pressure  distribution  with  constant 
arbitrary  upstream  pressure,  and  a  step  to  exactly  half  the 
upstream  pressure.  The  case  where  the  step  length  Is  2s,  l.e. 
a  full  wavelength,  corresponds  to  zero  upstream  pressure,  and 
constant  arbitrary  pressure  on  the  step. 

The  actual  free-surface  shape  generated  by  these  two 
pressures  Is  Illustrated  in  Figures  2.1  and  2.2  respectively. 
In  the  general  case  where  the  pressure  is  given  by  (2.8),  the 
free-surface  elevation  Is  given  by. (2.1)  as 

N 

n(x)  -  -P,+  £  (Pj+,-P j)K(x-Xj)  (2.12) 

Since  the  basic  kernel  function  K  is  easily  computed,  this 
provides  a  rapid  method  for  evaluating  linearised  solutions, 
not  only  for  the  present  case  where  P(x)  la  exactly  a  set  of  N 
steps,  with  N  small,  but  also  for  quite  general  P(x),  as  a 
numerical  procedure  in  which  (2.8)  is  an  approximation  for 
large  N.  Note,  however,  that  wherever  P(x)  possesses  a  finite 
jump  discontinuity,  n(x)  Is  finite  but  Its  derivative,  the 
free-surface  slope  n'(x)  is  infinite.  Specifically,  if  P(x) 
Jumps  by  AP  at  x  •  0,  then  near  x  ■  0  we  have 

n(x)  -  n(0)  -  ~  x  1  og  |  x  |  0  ( x )  (2.13) 

Figures  2.1  and  2.2  show  this  phenomenon  as  local 1 y- ver 1 1 c a  1 
free  surfaces  at  the  points  where  the  pressure  jumps. 

If  we  now  turn  to  the  case  N  ■  3,  since  there  are  now  two 
more  parameters  than  equations,  we  may  expect  to  obtain  a 
2-parameter  continuous  family  of  wave-less  solutions. 
Equations  (2.10)  now  state 


Linearized  free  surface  n(x)  generated  by  wave-less  pressures 
P(x)  that  vary  only  in  -1  <  a  <  0.  The  solid  curve  has  3  steps 
in  pressure,  the  dashed  curve  has  a  continuous  sinusoidally-varying 
pressure. 


P(x)  -  l-H(x)  -®* (x>  (2.16) 

where  H(x)  le  Che  unit  itep  function.  Thus,  fro*  (2.1),  we 

have 

n(x)  -  -1  -  K(x)  -  K" ( x ) 

-  -1  -(2co*  x-1)  -  (-2coax) 

•  0. 

Of  couree,  thie  Halting  eolutlon  le  very  singular  at  the 
origin,  and  for  exaaple  the  free-surface  elevation  becoaee 
unbounded  like  1/x,  but  It  Indicates  the  possibility  of 
wave-less  pressure  distributions  whose  domain  of  variation  is 
of  vanlshlngly-saall  length. 

The  logarithmic  singularity  (2.13)  in  the  free-surface 
slope  caused  by  Jumps  in  pressure  is  barely  acceptable  In  the 
linearised  theory,  and  becomes  (as  we  shall  see)  quite 
unacceptable  If  the  analysis  Is  extended  to  Include 
non-linearity.  It  Is  therefore  desirable  to  consider  families 
of  continuous  pressure  distributions.  The  simplest  such 
example  Is  a  linear  "pressure  ramp” 

Pi  ,  —  <  x  <  x,, 

P(x)  »  P,x/x,  ,  x,  <  x  <  0 

,0  x  >  0  ,  (2.1®) 

l/Lulv 

jfeatisf les(2.7)  If  x(  ■  -2v,  -4s,  -6s, .  Figure  2.4  shows 

the  free-surface  shape  at  x(  •  -2s.  This  curve  was  In  fact 
computed  using  (2.12)  with  N  ■  40,  which  gives  3-4  figure 
accuracy , 


(2.17) 


as  x  - 


A  useful  one-paraaeter  faally  of  continuous  pressures  is 


•Ivan  by 

P(x) 


Pj»ln 


ltx 


-•  <  x  <  xt 
x,  <  x  <  0 
x  >  0 


which  reduces  to  (2. IS)  if  P2  »  0,  but  Is  wave-less  for  mil 
x,  <  0  If 


(2.19) 


Figure  2.3  shows  (dashed  curve)  results  at  x,  •  -1.  The  solid 
and  dashed  curves  of  this  figure  contrast  discontinuous  and 
continuous  wave-less  pressures  for  the  sane  value  of  xt,  and 
clearly  the  saoother  pressure  gives  a  saoother  but  otherwise 
quite  similar  free  surface. 


If  one  demands  an  even  greater  degree  of  smoothness,  then 
for  example 


Kx)  • 


p‘*in,[£h-]  •  *«  <  *  < 


— •  <  x  <  x, 
0 

x  >  0 


(2.20) 


la  a  continuous  prssaure  distribution  with  continuous 
derivative,  and  Is  wave-less  If  x ,  •  -3v,  -Sr,  .... 

2.3  shows  results  at  xt  ■  ,jf t 


Figure 


Linearized  free  surface  n(x)  generated  by  a  wave-less  pressure 
P(x)  that  is  continuous  and  has  continuous  first  derivative, 
varying  sinusoidally  on  the  interval  On  <  *  <  0. 


3.  Soa a  propartita  ana  interpretations , 


The  linearised  free-eurface  flow  produced  by  e  weve-lese 
pressure  distribution  has  sose  Interesting  special  properties. 
In  particular.  symmetrical  wave-less  pressures  produce 
symmetrical  flows.  In  the  following  sense.  Suppose  x  ■  xs  Is  a 
point  of  syaaetry  of  a  wave-less  pressure  P(x),  such  that 


P(x,  *  X)  -  P(x,  -  X) 


(3.1) 


for  all  X.  Then  x  ■  x,  la  also  a  point  of  syaaetry  of  the 
free-surface  elevation  n(x)  produced  by  this  pressure,  l.e. 


n(x,  ♦  X)  -  n(x,  -  X) 


(3.2) 


Plgure  2.2  Illustrates  this  property,  with  xa  -  -s  as  the 
point  of  syaaetry. 


A  similar  result  applies  to  what  we  alght  call  "anti¬ 
symmetric"  wave-less  pressure  distributions,  but  Is  most 
directly  expressed  In  terms  of  derivatives.  That  Is,  If  x  •  x, 
Is  a  point  of  symaetry  of  the  pressure  gradient,  such  that 


P»(x,  ♦  X)  -  P'(x,  -  X) 


(3.3) 


for  all  X,  then  It  Is  also  a  point  of  syaaetry  of  the  fr« 


surface  slope,  l.e. 


n'(x,  ♦  X)  ■  n'(x,  -  x) 


(3. A) 


■ 

v  ’jlV^  v  ‘A 
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Upon  Integration,  w e  can  re-phrase  thla  intl-ayaattry  property 
to  atata  that 


P(x,  ♦  X)  -  P(x.)  -  P(x«)  -  P(x.  -  X)  (3.5) 

lnpllea  that 

n(x,  X)  -  n(xa)  ■  n(xa)  -  n(x,-  X)  (3.6) 

Figures  2.1,  2.3,  2.4  and  2.3  all  llluatrate  thla  property,  for 
varloua  valuea  of  xs. 

It  la  laiportant  to  eaphaalie  that  theae  propertlee  hold 
only  for  wave-leaa  preaaure  dlatrlbutlona.  It  la  obvioua  that 
aolutiona  with  wavea  cannot  obey  auch  ayaaetry  condition* , 
alnce  the  radiation  condition  denanda  that  the  wavea  occur  only 
at  x  "  +« ,  and  not  at  x  •  The  aynnetry  reaulta  can  be 

proved  by  noting  that  K(x)  in  (2.2)  can  be  epllt  Into  even  and 
odd  parta,  but  that  the  even  part  (coa  x  -  1/2)  contrlbutea  at 


■oat  a  conatant  to  n,  for  wave-leaa  preaaurea  P(x).  On  the 
other  hand,  the  odd  part  of  K(x)  neaearlly  generate*  (for  any 


P(x))  aolutiona  n(x)  aatlafylng  the  aynnetry  condition*. 

Another  alnoat  aelf-evldent  property  of  wave-leaa 
preaaurea  la  that  they  are  rerer»lbJe.  That  la.  If  P(x) 
goneratea  n(x),  then  P(-x)  generates  n(-x).  In  effect,  we  have 
alaply  reversed  the  direction  of  the  free  strean  at  infinity, 
and,  with  It,  the  whole  flow. 

Again,  thla  property  can  be  proved  fornally  for  the 
linearised  solution  by  separating  K(x)  Into  even  and  odd  parta. 


It  la  obviously  /else  for  wave-generst inn  pressures,  since  the 
waves  due  to  P(x)  must  lie  downstream,  which  means  that  those 


duo  to  P(-x)  would  lie  upstream  If  reversibility  held.  In 
fact,  this  reversibility  property  holds  not  only  for  weak 
pressures  P(x)  for  which  linearisation  Is  justified,  but  also 
in  the  fully-non-llnear  case. 

Once  a  flow  due  to  a  given  pressure  distribution  P(x)  has 
boon  computed.  It  may  be  re-interpreted  as  that  produced  In  the 
prasencs  of  a  certain  body.  That  Is,  the  computed  curve 
y  "  n(x)  is  a  streamline,  and  part  or  all  of  It  may  be  replaced 
by  a  fixed  Impermeable  boundary.  In  particular,  suppose 
F(x)  ■  0  for  x  <  xl  and  for  x  >  xg,  but  P(x)  a  0  for 
Xf,  <  x  <  xg.  Then  the  flow  generated  by  this  P(x)  Is  the  same 
as  that  past  a  body  with  the  equation  y  »  n(x),  x^  <  x  <  xg  , 
with  a  from  surface  at  atmospheric  pressure  In  x  <  xt,  and 
x  >  xg.  This  Interpretation ' holds  whether  or  not  the  pressure 
P(x)  Is  wave-less. 

In  this  way  we  can  provide  an  -invmrma  solution  to  problems 
Involving  ship-llke  bodies.  Examplss  such  as  that  in  Fig.  2.2 
whore  P(x)  la  non-sero  only  over  a  finite  segment  appear  at 
first  glance  to  be  the  most  relevant.  However,  pressures  P(x) 
that  are  non-sero  on  a  seml-lnflnite  segment,  such  as  that  used 
for  Fig. 2.1,  generate  flows  that  can  be  Interpreted  directly  as 
thoea  near  the  stern  of  a  ship.  In  general,  such  stern  flows 
generate  downstream  waves.  However,  In  the  special  case  where 
P(x)  Is  wave-less,  the  above  reversibility  property  applies, 
and  allows  an  additional  Interpretation  as  a  flow  near  the  bow 
of  a  ship. 


-20- 


S  '-'  ■.*  •*,  .  .  .  .  .  , 


.N 


'irhlftfi  "  nv^'~‘  ^"i  i~*~  v  ^ 


Of  course,  no  guarantee  la  available  that  the  "body"  shapes 
y  ■  n ( x )  so  generated  will  be  particularly  shlp-llke,  and 
finite-length  bodies.  Ilka  that  of  Fig,  2,2  tend  to  be  very  far 
froa  shlp-llke.  On  the  other  hand,  the  seal-infinite  body 
shown  In  (for  example)  Fig.  2.S  Is  not  unlike  a  bulbous  bow, 
and  tha  Inverse  technique  appears  to  be  quite  promising  If 
confined  to  local  flow  properties  at  the  bow  or  stern. 

In  asking  this  Inverse  Interpretation  of  the  linearised 
results  of  Section  2,  we  aust  not  forget  that  these  are 
approxlaat  Ions  ,  and  are  formally  only  valid  when  P(x)  Is  so 
weak  that  the  free-aurface  slope  n'(x)  Is  everywhere  small.  If 
y  ■  n(x)  la  now  Interpreted  as  the  surface  equation  of  a 
shlp-llke  body,  this  places  severs  restrictions  on  the 
applicability  of  the  results.  Fortunately,  tha  Inverse 
interpretation  Is  not  restricted  to  linearised  solutions,  and 
we  now  turn  our  attention  to  the  full  non-linear  problem,  whose 
solutions  allow  finite  n'(x),  and  hence  finite  slopes  In  the 
Inversely-generated  body  surface. 


-  '  < 


4. 


Less-Weak 


Pressures . 


Ir 


If  e  It  4  aetiuri  of  the  ill*  of  tho  pressure  P(x),  then 
the  results  of  Section  2  apply  In  the  Halt  as  a  -  0 ,  and 
provide  (a)  estimates  of  the  free-surface  shape,  l.e.  neglect 
*11  teras  of  0(ez)  (nd  saaller.  It  Is  nov  natural  to  seek  an 
expansion  In  powers  of  a,  commencing  with  the  0(e)  teras 
already  worked  out. 

However,  for  a  nuaber  of  reasons.  It  Is  Inconvenient 
(soaetlaes  impossible)  to  carry  out  this  analysis  for  a 
pressure  P  that  la  given  as  a  function  of  the  space  variable 
x.  Indeed,  if  the  pressure  is  large  enough,  there  aay  exist  no 
such  single-valued  functional  representation,  since  there  may 
be  aore  than  one  point  on  the  free  surface  corresponding  to  a 
given  value  of  x.  It  la  preferable  to  represent  the  pressure 
P  as  a  function  of  a  variable  that  changes  In  a  strictly 
aonotone-lncreaslng  Banner  along  the  free  surface. 

Such  a  variable  Is  the  velocity  potential  o.  That  Is,  we 
now  suppose  that  the  free-surface  pressure  P(e)  is  a  given 
function  of  e.  In  the  present  non-dlaens lonal 1  sat  Ion ,  the 
undisturbed  stream  corresponds  to  e  ■  x,  so  that  this 
representation  le  identical,  with  P(e)  e  P(x),  to  that  used  in 
Section  2,  to  the  order  of  approximation  being  used  in  that 
section.  All  results  of  Section  2  can  now  be  carried  over, 
with  x  replaced  by  e. 


We  now 
saae  time. 


consider  0<s*)  and  higher  perturbations.  At  the 
It  Is  convenient  to  suppose  that  the  pressure  P 


Hg 

& 

1*9 
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Itself  can  be  modified  as  we  expand  Co  such  higher  orders. 


P(e)  •  eP, (a)  +«,p,(»)  +  «'Pj(») 


<4.  l> 


Since  e  was  scaled  to  unity  In  Section  2t  In  fact  the  results 
of  that  section  correspond  to  Pt(«)  ■  P(x). 


If  s  •  xely,  and  f  ■  e  ♦  !♦  la  the  complex  velocity 
potential,  then  we  consider  s  “  *(f)  and  expand 


■  <f>  ■  f  ♦  ai,(f )  +  a 2  x ,  ( f )  +  e*tl(f) 


(4.2) 


The  free  surface  corresponds  te  t  •  0,  and  the  (scaled)  dynamic 
free-eurface  condition  Is  (with  g  ■  1 ,  U  •  1,  p  ■  1) 


P  ♦  | l • (f ) | “*/2  ♦  y  •  1/2 


(4.3) 


Upon  substitution  of  (4.1)  and  (4.2)  and  collecting  powers  of 


a,  we  find  that  for  all  Jal,2,3. 


Pj+pj  -  *j'  ♦  yj  ■  0 


(4.4) 


where  all  quantities  are  evaluated  on  e  •  0,  the  prime  denotes 


3/a*,  and 


Pt  •  0 


3  <  t  l  •  t 

Pt  •  — *  \  m  — y  | 

2  2 


ps  -  3x ,  *,  -  y,  y t  +  y»  -  2x, 

p 4  •  *3  *  yi  y»  ♦  ***  -  ?y* 


(4.5) 


(4.6) 


(4.7) 


,  I  ,  »  •  I  •  -  I  •  J  . 

-  6x ,  'x,  ♦  4  x  j  y,  y,  ♦  2x,  y,  *-  lx, 

-  I*.**  -  }y,*‘ 


*  y , '  * 

(4.8) 


etc.  Thus,  st  all  orders  of  ipproxlastion  J,  ws  have  the  sasie 
type  of  boundary  condition  (4.4)  the  only  complication  at 
hither  order  being  the  aanner  In  which  the  "equivalent  linear" 
pressure  Pj  Is  determined  from  earlier  terns  In  the 
expansion. 


For  J  ■  1,  (4.4)  and  (4.5)  are  equivalent  to  the  usual 
linearised  free-surface  boundary  condition,  and  all  results  of 
Section  2  are  available  to  yield  solutions  of  this  probles. 
Note  that  In  the  present  notation. 


yt(«)  "  n(e) 


(4.9) 


In  fact,  we  need  no  wore  Information  from  Section  2  than  n(x) 
(with  x  replaced  by  e),  since  (4.4)  and  (4.9)  imply 


*»'(•)  -  Pt(#)  ♦  n(e> 


(4.10) 


If  we  turn  to  the  2nd-order  terns  J  ■  2,  It  Is  lanedlately 
apparent  that  step-function  discontinuities  In  ths  Input  P,(e) 
are  unacceptable ,  since  they  lead  to  (c.f.(2.13)) 
logarlthmlcally-lnf lnlte  free-surface  slope, and  via  (4.6)  to 
"log  squared"  singularities  In  the  2nd-order  pressure  p2  ,  and 
hence  to  an  unbounded  2nd-order  free-surface  displacement  yt  . 
This  problen  compounds  Itself  at  higher  order,  and  one  must 
suspect  that  no  solution  exists  if  P(e)  possesses  a  step 
discontinuity. 
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Thlo  la  alao  a  concluaion  one  la  drawn  to  by  direct 
conelderatlon  of  the  non-linear  boundary  condition  (4.3), 
lrreepectlve  of  any  form  of  aaynptotlc  expanelon  auch  aa  (4.2), 
and  appeara  to  be  In  conflict  with  the  fact  that  Schwartz  [7] 
clalaa  to  have  found  a  non-linear  nuaierlcal  aolutlon  for  auch  a 
preaaure  dlatr Ibutlon.  However,  It  la  probable  that  Schwartz'a 
nuaierlcal  aolutlon  la  not  aufflclently  accurate  near  the 
alngular  point*  to  display  thla  difficulty. 


In  any  caae,  we  chooae  to  require  that  Pt(o)  be  a 
continuous  function,  typically  (2. IS),  (2.19)  etc.  For 
exaaple.  If  P,  la  contlnuoua  but  Ita  derivative  P,'  poaaeaaea  a 
atep-f unct Ion  discontinuity,  then  y, 1  -  nt’  la  bounded  and 
contlnuoua,  and  hence  ao  la  the  aecond-order  preaaure  p2.  We 
therefore  expect  no  trouble  at  any  order  of  approxlastlon,  and 
by  laiplicatlon,  expect  that  the  full  non-llnear  problea  doea 
have  a  aolutlon. 


The  ayaiaietry  propertlea  of  the  linear  wave-leaa  aolutlona 
aa  dlscuaaed  In  Section  3,  now  play  an  Important  role  In  the 
non-llnear  theory.  Suppoae  for  exaaiple  that  Pt(e)  la  wave- 
leaa  and  ayaaetrlc  about  aone  point  »  •  *(.  Then,  ao  we  have 
aeon,  ao  la  n(e)  ■  y,(e),  and  thus  ao  la  xt'(e).  Hence 
Pz  -  Pj  la  aynnetrlc,  and  a  ayaaetrlc  aecond-order  aolutlon 
exlata  If  the  aecond-order  Input  P(  la  taken  aa  aynnetrlc.  A 
alnllar  concluaion  appllea  at  all  higher  orders,  and  we  nay 
therefore  expect  to  generate  a  fully  non-llnear  aynnetrlc 
aolutlon.  Note  that  In  auch  a  caae  P,  la  not  arbitrary  and 
cannot  be  aet  to  zero;  It  la  needed  In  order  to  cancel  the 
wavea  at  aecond  order. 


The  eti*  of  "antl-ayaaetrlc”  preaaurea  la  aoaewhat  ter  a 
Interact  t  ln| .  Now  if  cha  praaaura  gradient  P,'  la  ejvaaatrlc, 
than  ao  la  tha  fraa-aurfaca  alopa  n'  ■  yt ' .  Thua  (3.6)  lapllea 
that  p}  -  Pt  la  ayaaatrle,  and  wa  can  (lndaad  tuat)  uaa  a 
ayaaatrle  2nd-order  Input  Ft.  Thua,  at  2nd-ordar,  tha  antl- 
ayiaatry  proparty  of  tha  Input  praaaura  la  nacaaaarlly 
daatroyad,  and  It  appaara  that  thara  axlata  no  viva-laaa 
non-llnaar  ant 1-eyaae t r 1 c  praaaura  d 1  a t r lbut Ion . 

For  axaapla,  auppoaa 

{1,  -•  <  a  <  -2a 

-♦/2a,  -2a  <  e  <  0 

0,  a  >  0  (4.11) 

Thla  la  Juat  tha  "raap"  praaaura  of  (2.18),  expreaaed  In  taraa 


and 

•valuated 

In 

the  wava-laaa 

C  ASA 

x,  “  -2«r.  Note  that 

thla 

pressure 

Is 

antl-ayaaatrlc 

About 

a  ■  -a.  Thla  praaaura  can  ba  uaad  to  coaputa  x,(e)  and  yt(e), 
tha  lattar  bains  aa  plotted  for  n(x)  In  Figure  2.4.  Than  wa 
can  eaelly  coaputa  pt(e)  froa  (4.6),  and  thla  la  given  In 
Figure  4.1,  which  varlfiaa  that  Pt(e)  la  ayaaatrle  about 
a  w  -a.  Tha  coaplate  aecond-order  problaa  now  haa  bean  reduced 
to  a  linear  problaa.  In  which  tha  free  aurfaca  la  dlaturbad  by 
a  coablnatlon  of  tha  input  aacond-ordar  praaaura  Pf(e),  (which 
aay  or  way  not  ba  praaent)  and  tha  coaputad  praaaura  p«(a)  aa 
dlaplayad  In  Figure  4.1,  which  auat  alwaya  ba  praaent. 

If  thara  la  no  Input  aacond-ordar  praaaura,  l.a.  If  tha 
praaaura  la  given  by  a  alapla  rasp  to  at  laaat  0<a')  accuracy, 
than  wavaa  auat  ba  generated.  That  la,  by  ltaelf,  tha  praaaura 
Pi(e)  of  Figure  4.1  deer  generate  wave*.  If  wa  deaand  a 


Tit  solid  curve  is  the  secona-order  pressure  p  ($)  corresponding 
eo  sn  input  first-order  pressure  equivalent  to  the  linear  ramp  of 
Fig.  2.4.  The  dashed  curve  is  a  wave-less  analytic  approximation 
to  the  solid  curve  for  -2 it  <  4  <  0,  given  by  equation  (4.15). 


vtvflili  solution  at  aacond  order,  we  oust  Incorporate  some 


Input  preaaure  F,(«)  at  second  order,  of  such  a  nature  and 
Magnitude  as  to  cancel  the  waves  generated  by  p,(e)>  The 
oymmatry  principle  desands  that  this  input  pressure  must  be 
symmatrlc  about  o  “  -»,  and  there  are  many  choices  for  P,(e) 
that  are  satlfactory. 

A  rather  attractive  possibility  Is  to  choose 

P*(a)  a  -Pi (a)  (4.12) 

in  which  case  there  Is  no  pressure  distribution  at  all  at 
2nd-order,  and  hence  no  second-order  flow.  That  Is,  without 
further  ado,  we  can  assert  that  an  Input  pressure 


*<•>  •  eP,(a>  -  e'p,(e)  (4.11) 

•here  Pt(a)  Is  given  by  (4.11)  and  Pt(e)  by  Figure  4.1  produces 
a  wave-lass  flow  field. 

s  •  ««,(!)  ♦  0(a  *)  (4.1*) 

One  difficulty  of  the  assumption  (4.12)  is  that  since 
P«(e)  S  0  In  a  >  0,  the  same  applies  to  Ps(e).  That  Is,  there 
la  no  true  free  surface  on  which  -p  »  0,  the  Input  pressure 
being  non-aero  everywhere. 


Alternatively,  a  good  choice  for  P,(a)  is  one  that,  while 
remaining  taro  for  a  >  0  and  generating  exact  ly  equal  and 
opposite  waves  to  pt ,  also  cancels  p>(e)  for  -2v  <  a  <  0,  to  as 


grace  an  extant  aa  possible.  One  auch  choice,  found  by  trial 
and  error, 1  a 

0.118  a  In  3o/2  4  0.09  eln*e,  -2a  <  •  <0. 

-P*  <♦)  - 

0  ,  otherwlee 

(4.15) 

which  la  ahown  (daahed)  In  Figure  4.1.  Kota  that  the  flrat 
ten  of  (4.15)  generatea  the  wavea  to  cancel  thoae  of  p2;  the 
aecond  term  la  itaelf  wava-laaa,  and  the  coefficient  0.09  la 
alaply  choaen  to  natch  -Pi  to  pt  aa  eloeely  ae  poaalbla.  It  le 
now  quite  a t r eight - f orward  to  compute  the  free-aurfaca 
elevation  y2  produced  by  the  ever ywhere-anal 1  residual  pressure 
P 2+p g ,  and  this  Is  shown  In  Figure  4.2.  Note  the  snail 
absolute  size  of  v ,,  which  suggests  that  these  results  may 
retain  validity  up  to  quite  large  values  of  a. 

The  actual  free-aurface  shape  In  physical-space 
coordinates  (x,y),  is  given  by 

x  •  a  *  ex, (e)  +  0(e  *  ) 

y  “  ay, (a)  +  e*y*(e)  *  0(a*>.  (4.16) 

To  obtain  a  result  for  y  ■  y(x),  consistently  with  error  0(e*)» 
we  can  treat  (4.16)  as  a  pair  of  paranatrlc  equations,  with  a 
as  a  paraneter.  The  small  magnitude  of  y*(a)  In  the  example  of 
Figure  4.2  means  that  (for  moderate  e),  the  y-coordlnate  of  the 
free  surface  is  little  changed  by  second-order  effects  In  this 
example.  However,  the  0(e)  correction  to  the  x-coordinate  is 
not  so  small,  and  represents  the  most  significant  non-linear 
modification.  We  return  to  thla  question  in  the  following 


sect  Ion 


3.  Strong  fr«*«ur«« 

An  *x«et  inverts  method  of  solution  of  ths  present  close  of 
problea  is  always  ovoileble.  That  is,  given  eny  function  z(f) 
that  la  analytic  in  the  lower  half  f-plane,  we  can  use  (3.3) 
directly  to  coaputa  the  free-surface  preeeure  P(e).  One  way  to 
do  this  le  to  coablne  a  unlfora  etreaa  with  singularities  above 
the  free  surface. 

For  exaaple,  suppose 

»<f>  •  t  *  12  log(f-Ai)  ♦  0.1) 

If  linearisation  (i.e.  fat)  were  legltiaate,  (3.1)  could  be 
Interpreted  as  a  etreaa  plus  a  oink  at  (0,4)  and  a  dipolo  at 
(-4,12).  However,  the  else  of  the  coefficients  in  (3.1)  makes 
linearisation  of  doubtful  value,  even  for  qualitative  purposes. 
Figure  3.1  shows  the  streaallne  x+iy  -  z(e+10)  coaputed  froa 
(3.1).  Since  we  have  scaled  g  •  U  ■  1,  and  the  "draft"  of  the 
shape  shown  in  Figure  3.1  is  12*,  the  draft-based  Froude  nuaber 
for  this  flow  la  (12*)”*^*  ■  .143.  The  shape  shown  is  quite 
like  a  typical  bow. 

The  preeeure  F  (obtainable  froa  (3.3))  is  alaost 
hydrostatic,  and  is  significant  only  for  y  <  0.  Indeed,  if 

y  >  0,  then  P/P,  <  .008,  where  P,  is  the  pressure  at  x  - 
The  coefficients  and  singularity  location*  were  in  fact  chosen 
by  trial  and  error  so  that  this  would  be  so.  Many  other 

choices  of  the  parameters  in  (3.1)  led  to  a  similar  property, 
and  shapes  having  snail  pressure  for  sufficiently  lerge  e  or  x 
can  be  constructed  at  any  draft-based  Froude  nuaber.  However, 
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this  procedure  seems  hard  Co  carry  through  in  any  sy 


ic 


manner,  and  we  turn  to  an  alternative  procedure  based  on  the 
linearized  solution. 

Flows  satisfying  the  exact  free-surface  condition  (3.3)  can 
be  constructed  in  a  systematic  manner  by  making  use  of  the 
results  already  obtained  in  the  previous  sections.  For 
example,  suppose  we  set 

z  -  f  e  z , <  f )  (5.2) 

(exactly),  where  z,(f)  is  the  linear  solution  corresponding  to 
some  input  "seed"  pressure  Pt(e).  Then  the  flow  (3.2)  can  be 
considered  to  have  been  generated  by ■ the  free-surface  pressure 
p(e)  given  exactly  by 

p(e)  -  X  -  e  y , ( e )  -  5  1(1  ♦  ex, * (e) ) *♦  (ey, • (e))*  ]** 

*  (5.3) 

Note  that  in  order  to  generate  the  free-surface  shape  and 
pressure  P,  we  need  only  compute  the  linearised  free-surface 
elevation  n(x)  with  e«l,  since  then  jr,(«)  ■  n(e),  and 
Xt'(e)  •  Pi(e)+n(e).  We  can  then  compute  very  rapidly  a 
complete  family  of  exact  non-linear  solutions,  by  varying  the 
amplitude  parameter  e.  The  free-surface  shape  is  given 

peremetr ically  by 


S3 

& 

Jv 


3 


*•••*■  eJ*-IP,(e)+n(e)  J  de 

y  •  en(e)  (5.4) 

As  e  varies,  the  y-coordinate  simply  scales  in  proportion,  but 
the  actual  shape  of  the  free-surface  distorts,  because  of  the 
non-proportional  dependence  of  x  on  e . 


$ 
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Figure  3.2  shows  free-eurface  shapes  eosputed  in  this  way 
for  various  values  of  a,  with  Che  seed  pressure  Ft(a)  given  by 
the  slaiple  ramp  (4.11).  These  curves  are  drawn  true  to  scale, 
and  their  non-linearity  Is  confined  by  the  large  slopes. 
However,  the  curve  shown  for  eel  does  not  display  large 
slopes,  and  hence  it  la  almost  the  sane  as  Figure  2.4,  since 
linearisation  is  not  too  much  In  error.  As  a  Increases,  the 
"downward  bulb"  at  a  ■  -3.3  becomes  sharper,  whereas  the  "upwrd 
bulb"  at  e  ■  -1  (which  Is  in  Figure  2.4  exactly  a  reversed 
Image  of  the  downward  bulb)  is  smoothed  out. 


Eventually,  at  about  e  «  3.1,  the  downward  bulb  sharpens 
Into  a  cusp,  and  for  a  values  greater  than  this,  non-physical 
looped  shapes  are  generated.  The  -distorted  variation  in  the 
x-coordlnate  now  allows  a  non-single-valued  y  •  y(x)  free 
surface,  containing  a  point  where  the  freo  surface  la  vortical. 


and  this  actually  happens  for  a  values  near  to  3  in  this 
example . 


For  each  separate  value  of  e,  wo  must  compute  the  actual 


free-surface  pressure  P(e)  by  (3.3). 


Figure  3.2  corresponds  to  a  different  P(e).  As  the  downward 
bulb  becomes  sharper,  the  fluid  velocity  lncreaeee  near  the 
corner,  and  hence  the  pressure  F(a)  decreases,  until  at  the 
value  of  e  where  the  body  Is  cusped,  F  -  -■  at  the  cusp. 
Figure  3.3  (solid  curves)  repeats  the  froe-surface  shape  of 
Fig.  3.2  at  e  •  2,  and  shows  the  corresponding  pressure  F, 
plotted  as  a  function  of  x. 
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A  Mt  Of  oxact  non-linear  streanlines  </  ■  0  for  flows  seeded  by 
the  linear  raap  pressure  of  Fig.  2.4,  for  various  values  of  a 
paraaeter  c  Measuring  the  Magnitude  of  the  input  pressure. 
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In  fact  It  Is  Interesting  to  note  that  the  present  exact 
procedure  generates  alaost  the  saae  results  as  the  2nd-order 
procedure  (4.16).  The  only  difference  between  (4.16)  and  (5.4) 
Is  the  *  term  e*y,(«)  In  the  former,  and  we  have  Indicated  In 
Section  4  that  this  tera  can  be  aade  quite  small.  That  Is,  the 
second-order  accurate  free  surface  computed  at  a  •  2  using 
(4.16),  with  Pj  given  by  (4.11)  and  P*  by  (4. IS),  is  alaost 
indistinguishable  from  that  shown  in  Figure  S.3.  However,  the 
pressure  P(o)  •  eP,  +  e*P,  now  vanishes  exactly  for  all  0  >  0, 
so  that  there  is  a  true  free-surface  present. 

The  present  exact  Inverse  nethod  can  be  extended  to 
account  of  both  1st  and  2nd-order  seed  preesure. 
suppose  Instead  of  (S.2)  we  set 

*  ■  f  *■  ez,<f )  +  «*t,(f)  (5.5) 

where  s,(f)  and  xt(f)  are  1st  and  2nd-order  solutions  generated 
by  P  ,  (♦)  and  Px(e)  respectively.  Then  this  Is  also  the  exact 
solution  generated  by  the  pressure 

P(b)  •  -ey,(e)  -  e*yx(e) 

-  jBl+ex, '(♦)+*,x,,(#))»*(«yi'(*)+«,yi '<•))*  ]“' 

(5.6) 

Figure  5.3  also  shows  (dashed)  computations  using  (5.5),  (5.6), 
for  e  •  2,  with  P,  given  by  (4.11)  and  Px  by  (4.13).  The 
dashed  end  solid  curves  are  generally  close,  although  the 
dashed  downward  bulb  Is  significantly  sharper,  and  hence  the 


take 

That  Is, 


However,  the  most  important  change  is  not  eesy  to  see  on 
the  scale  of  the  figure,  and  this  is  the  fact  chat  there  is  e 
range  (x  >  -1.4)  over  which  the  pressure  P  computed  from  (5.4), 
satisfies  |  P  |  <  0.01.  In  the  sane  range  x  >  -1.4,  the  solid 
curve  still  has  only  |P|  <  0.06.  That  is,  this  portion  of  Che 
dashed  streamline  is  much  more  nearly  "free"  then  the  solid 
one.  Clearly  one  could  continue  this  process  to  even  higher 
orders,  generating  each  time  an  exact  solution,  with  the 
property  that,  as  the  order  of  the  process  increases,  P  -  0 
for  e  >  0. 


6.  Conclusion 

We  have  demonstrated  in  this  paper,  both  bp  linear  and 
non-linear  methods,  some  families  of  free-eurface  pressure 
distributions  that  do  not  generate  waves.  Since  these 
wave-leas  flows  can  be  reversed  in  direction,  they  also 
generate  candidate  shapes  for  splash-less  bows. 
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